INPUT 3D 
COORDINATES 
OF ALL ATOMS IN A 
SELECTED MOLECULE 



FOR EACH ATOM, 
DETERMINE THE 
PREDICTED GIBBS 
FRE ENERGY OF 
BINDING TO AN IDEAL 
UGAND 



GENERATE 3D MODEL OF 
MOLECULE & MAP 
CORRESPONDING 
PREDICTED GIBBS FREE 
ENERGY OF BINDING TO 
EACH ATOM 



OUTPUT THE 3D 
MODEL AS A 
PREDICTION MODEL OF 
BINDING TARGETS 



FIG.1 



INPUT ATOMIC COORDINATES OF ALL ATOMS IN 
TARGET MOLECULE 
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CALCULATE BINDING POTENTIAL OF EACH ATOM IN 
TARGET MOLECULE 
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IDENTIFY REGIONS WITH HIGHEST BINDING POTENTIAL 



IDENTIFY POTENTIAL BINDING SITES IN TERMS OF 
BINDING POTENTIAL 



CLASSIFY POTENTIAL BINDING SITES ACCORDING TO 
CHEMICAL NATURE (E.G., DEGREE OF HYDROPHOBICS, 
CHARGE, ETC.), SURFACE AREA, LOCATION IN CAVITY, ETC. 



SELECT POTENTIAL BINDING SITE FOR UGAND DESIGN 

ACCORDING TO USER CRITERIA, E.G.: 

-TYPE OF UGAND (E.G., PEPTIDE, ORGANIC MOLECULE) 

-BINDING AFFINITY 

-SIZE 

-GEOMETRY 
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FIG.3 





CHARACTERIZATION OF BINDING SITE Wfffl RESPECT TO: 
-BINDING POTENTIAL 
-CHEMICAL NATURE 
-GEOMETRY 
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SELECTION OF APPROPRIATE LIBRARY OF SEED DIPEPT1DES: 
-400 DEPEPTIDES OF NATURAL AMINO ACIDS 
-SPECIAL LIBRARIES OF NON-NATIVE OR CHEMICALLY MODIFIED 
AMINO ACIDS 
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COARSE GRAINED SORTING OF SED DIPEPT1DES BASED UPON 
GEOMETRICAL AND CHEMICAL CHARACTERISTICS OF BINDING SITE 
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PLACE SEED DIPEPTIDE ON BINDING SITE 
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ENERGY MINIMIZE BOUND DIPEPTIDE CONFORMATION BY ROTATION 
AROUND BACKBONE AND SIDE CHAIN TORSIONAL ANGLES 










CALCULATE EXPECTED BINDING AFFINITY FOR DIPEPTIDE 
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SELECT BEST DIPEPTIDE SEEDS AS THOSE WITH THE HIGHEST 
BINDING AFFINfflES 
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BUILD LEAD PEPTIDE UGANDS STARTING FROM SELECTED SEED 
DIPEPTiDES BY ADDITION OF AMINO ACIDS AT EITHER END OF THE 
CHAIN (REPETITION PROCEDURE) 
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ENERGY MINIMIZE THE CONFORMATION OF THE LEAD PEPTIDE AFTER 
EACH AMINO ACID ADDITION 
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CALCULATE EXPECTED BINDING AFFINITY 
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CONTINUE THE ADDITION OF AMINO ACIDS UNTIL NO INCREASE IN 
BINDING AFFINITY IS OBSERVED 
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SELECT LEAD PEPTIDE WITH HIGHEST BINDING AFFINITY 
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FIG.4 





STRUCTURE-BASED THERMODYNAMIC ANALYSIS OF COMPLEX OF 
LEAD PEPTIDE WITH MACROMOLECUUE 
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CALCULATION OF CONTRIBUTIONS OF INDIVIDUAL PEPTIDE RESIDUES 
OF BINDING ENERGETICS. CHANGES IN LEAD PEPTIDE WITH THE 
GOAL OF IMPROVED BINDING AFFINITY ARE USUALLY MADE AT 
THOSE POSITIONS THAT CONTRIBUTE THE LEAST TO GIBBS ENERGY 
OF BINDING 
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MODIFICATION OF LEAD PEPTIDE IS ACCOMPLISHED BY A 
COMBINATION OF ONE OR MORE OF THE FOLLOWING OPERATIONS: 

1. ) MUTATION 

2. ) ADDITION OF AMINO ACIDS AT AMINO AND /OR CARBOXY 
TERMINUS OF LEAD PEPTIDE 

3. ) DELETION OF AMINO ACIDS AT EITHER END OF LEAD PEPTIDE 
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-AMINO ACIDS FOR MUTATION OR ADDfFlON ARE SELECTED FROM 
20 NATURAL AMINO ACIDS, NON-NATURAL AMINO ACIDS OR 
CHEMICALLY MODIFIED AMINO ACIDS 
-AMINO ACIDS ARE SYSTEMATICALLY CHOSEN AND INSERTED 
AT DESIRED POSITION FOR ADDITION OR REPLACEMENT 
-ATOMIC COORDINATES OF MODIFIED COMPLEX ARE GENERATED 


^8400 










FOR EACH CHANGE OR ADDITION THE CONFORMATION WITH THE 
LOWEST ENERGY IS IDENTIFIED BY SYSTEMATIC ROTATION AROUND 
BACKBONE AND SIDE CHAIN TORSIONAL ANGLES 
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EXPECTED BINDING CONSTANT IS CALCULATED 
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3D MODEL(l.E., ATOMIC COORDINATES) OF MODIFIED PEPTIDE 
BOUND TO MACROMOLECULE IS BUILT. 
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FIG.5 





CHARACTERIZATION OF BINDING SITE WITH RESPECT TO: 
-BINDING POTENTIAL 
-CHEMICAL NATURE 
-GEOMETRY 
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BUILD IDEAL UGAND BY PLACING ATOMS THAT ARE ENERGETICALLY 
COMPLEMENTARY TO THOSE IN THE BINDING SITE AT OPTIMAL 
VAN DER WAALS DISTANCES 
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SELECT LEAD COMPOUND USING GEOMETRICAL ARRANGEMENT OF 
ATOMS IN IDEAL UGAND AS A TEMPLATE. SELECTION IS MADE FROM 
DATABASE OF COMPOUNDS OR BY BUILDING MOLECULE WITH 
APPROPRKTE ATOMIC ARRANGEMENT. 
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I 1 PLACE LEAD COMPOUND ON BINDING SITE. 
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ENERGY MINIMIZE BOUND CONFORMATION BY SYSTEMATIC BOND 
ROTATION. 
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CALCULATE EXPECTED BINDING AFFlNflY FOR LEAD COMPOUND 
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SELECT BEST LEADS AS THOSE WITH THE HIGHEST BINDING 
AFFINITIES. 
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OPTIMIZE LEAD COMPOUND UGANDS BY ADDITION OR REPLACEMENT 
OF CHEMICAL GROUPS. 
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ENERGY MINIMIZE THE CONFORMATION OF THE LEAD COMPOUND 
AFTER EACH MODIFICATION. 
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CALCULATE EXPECTED BINDING AFFINITY. 
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CONTINUE MODIFICATION UNTIL NO INCREASE IN BINDING AFFINITY IS 
OBSERVED. 
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SELECT LEAD COMPOUND WITH HIGHEST BINDING AFFINITY. 
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FIG.6 



STRUCTURE-BASED THERMODYNAMIC ANALYSIS OF COMPLEX OF 
LEAD COMPOUND WITH MACROMOLECULE 



CALCULATION OF CONTRIBUTIONS OF INDIVIDUAL ATOMS TO BINDING 
ENERGETICS. CHANGES IN LEAD COMPOUND WITH THE GOAL OF 
IMPROVED BINDING AFFINITY ARE USUALLY MADE AT THOSE 
POSITIONS THAT CONTRIBUTE THE LEAST TO GIBBS ENERGY OF 
BINDING 
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MODIFICATION OF LEAD COMPOUND IS ACCOMPLISHED BY A 
COMBINATION OF ONE OR MORE OF THE FOLLOWING OPERATIONS: 

1. ) REPLACEMENT OF CHEMICAL GROUPS 

2. ) ADDITION OR DELETION OF CHEMICAL GROUPS 







-CHEMICAL GROUPS FOR REPLACEMENT OR ADDITION ARE 
SELECTED FROM TOOLBOX OF STANDARD ORGANIC GROUPS: 
E.G., METHYL, AMINO, HYDROXYS ETC. 
-AFTER CHEMICAL MODIFICATION OF LEAD COMPOUND, 
THE ATOMIC COORDINATES OF THE MODIFIED COMPLEX 
ARE GENERATED. 
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FOR EACH CHANGE OR ADDITION THE CONFORMATION WITH THE 
LOWEST ENERGY IS IDENTIFIED BY SYSTEMATIC BOND ROTATIONS 






EXPECTED BINDING CONSTANT IS CALCULATED 






3D MODEL (I.E., ATOMIC COORDINATES) OF MODIFIED COMPOUND 
BOUND TO MACROMOLECULE IS BUILT. 
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FIG.7 
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FIG.10B 
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FIG.10C 
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FIG.15B 
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FIG.19B 



